REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reDortinq  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  Instructions,  searching  existing  data  sources, 
aatherina  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
Mllectlon  of  Information  inciudinq  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports.  1215  Jefferson 
Davis  Highway  Suite  1204  Ariin^on  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington.  DC  20503. 

1  AGENCY  USE  ONLY  I  2.  REPORT  DATE  I  3.  REPORT  TYPE  AND  DATES  COVERED 

3/17/98 _ Final  9/1/93-8/31/97 _ 


4.  TITLE  AND  SUBTITLE 


First  Principles  Calculations  Using  Linear  Response 
Theory  for  Ferroelectrics 


5.  FUNDING  NUMBERS 

Contract  No.: 
K00014-93-1-1095 


6.  AUTHOR(S) 

Henry  Krakauer 


7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(ES) 

College  of  William  and  Mary 
Department  of  Physics 
P.  0.  Box  8795 

Williamsburg,  VA  23187-8795 

9.  SPONSORING  /  MONITORING  AGENCY  NAMES(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 


19980320  047 


a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


1  mwution  code 


BWiugtmOR  BTATE^eWT  g 

AjipioTod  ias  public  tsiMMi 
DjiUlbnttoe  _ 


13.  ABSTRACT /Max/mom  200  words/ 


First-principles  density- functional  methods  using  LAPW -and  LAPW  linear  response 
method  were  applied  to  study  the  incipient  ferroelectric  SrTi03.  Despite  similari¬ 
ties  to  BaTi03,  SrTi03  remains  paraelectric  down  to  zero  .temperature ,  transforming 
instead  to  an  antiferrodisto'rtive  structure  at  IO5  K.  As  a' function’ Of  wavevector , 
the  linear  response  calculations  show,  structural  instabilities  along  the  R-M-R  line 
in  the  Brillouin  zone,  consistent  with  this  antiferrodistortive  ground  state. 
Although  a  ferroelectric-type  instability  was  also  found  near  the  zone  center  as  in 
BaTi03,  the  phase  space  of  this  instability  is  greatly  reduced  by  comparison. 
Moreover,  we  find  that  the  antiferrodistortive. phase  is  stable,  marginally,  against 
ferroelectric  structural  distortions.  In  combination,  these  findings  help  explain 
the  absence  of  a  transition  to  a  ferroelectric  structure. 

DTIC  QUALITY  mSPECTED  6 


14.  SUBJECT  TERMS 

Ferroelectrics,  first-principles  calculations,  density 
functional  theory,  LAPW  method,  linear  response 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

unclassified 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

unclassified 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

unclassified 


20.  LIMITATION  OF  ABSTRACT 


standard  Form  298  (Rev.  2- 
Prescribed  by  ANISE  Sad  Z39-18 
298-102 


Final  Report  on  ONR  AASERT  Grant  N0001 4-93-1 -1095  —  PI:  Henry  Krakauer,  College  of  William  &  Mary  Page  2 

A.  Description  and  Motivation  of  the  Scientific  Research  Goals 

Although  it  possesses  the  same  high-temperature  cubic  perovskite  structure  as 
ferroelectric  materials  like  BaTiOs  and  KNb03,  SrTi03  exhibits  qualitatively  different 
temperature-dependent  behavior.  It  transforms  to  a  non-ferroelectric  phase  when  the  temperature 
is  lowered  below  Tc  =  105  K,  due  to  a  soft-mode  antiferrodistortive  (AFD)  instability  at  the  R- 
point,  the  Brillouin  zone-comer.  The  atomic  motions  in  this  AFD  mode  are  characterized  by 
rotations  of  the  TiOe  octahedra  about  a  [1  0  0]  axis.  Below  Tc,  these  octahedra  are  slightly 
rotated,  with  the  twist  angle  about  the  tetragonal  c-axis  being  the  order  parameter  of  the 
transition.  Upon  further  cooling,  the  dielectric  constant  rises  in  a  Curie-Weiss  manner  until  about 
T  =  35  K,  where  it  continues  to  increase,  but  in  a  non-singular  manner,  and  at  T  =  3  K  and  below, 
it  saturates  to  a  value  of  about  10“^,  but  a  further  stmctural  phase  transition  does  not  occur.  Due  to 
these  unusual  properties  at  low  temperature,  SrTiOs  has  been  subjected  to  extensive  studies  in 
the  past  few  decades.  In  recent  years,  first-principles  density  functional  theory  calculations  using 
the  local  density  approximation  (LDA)  have  achieved  remarkable  success  in  determining  the 
physical  properties  of  ferroelectrics,  including  the  equilibrium  atomic  configurations,  lattice 
vibrational  frequencies,  and  the  underlying  microscopic  mechanisms  behind  macroscopic 
behaviors.  For  instance,  our  recent  first-principles  linear  response  calculations  for  KNb03 
revealed  large  planar  regions  of  phonon  instabilities  in  the  Brillouin  zone.  These  correspond,  in 
real  space,  to  chains  along  the  [1  0  0]  directions  of  atoms  coherently  displaced  along  the  chain 
direction.  Subsequent  molecular-dynamics  simulations  of  KNb03  have  shown  that  these  chains 
are  dynamic  rather  than  static  and  that  they  are  preformed  even  in  the  high-temperature 
paraelectric  phase.  [2]  (Note:  references  in  square  brackets  refer  to  publications  in  Section  D.) 
Motivated  by  this  success,  we  undertook  a  first-principles  study  of  SrTiOs  in  both  cubic  and  AFD 
phases  to  examine  the  nature  of  its  stmctural  instabilities  to  explSin  the  its  distinctive  properties. 

B.  Significant  Results 

Comparing  with  experimental  phonon  frequencies  obtained  with  neutron  scattering  and 
infrared  spectroscopy  at  high  synunetry  points  in  the  Brillouin  zone,  the  theoretical  phonon 
frequencies  are  found  to  be  in  good  agreement.  [1, 3]  In  addition,  we  correctly  obtain  the  known 
instability  at  the  R-point  in  the  Brillouin  zone  (BZ),  which  is  responsible  for  the  AFD  phase 
transition  at  105  K.  Moreover,  we  found  that  the  AFD  instability  extends  along  the  entire  R-M-R 
line  in  the  BZ,  forming  a  one-dimensional  cylindrical  tube  in  the  BZ.  We  also  found  that  this 
instability  competes  with  a  ferroelectric  instability  centered  at  q  =  0,  whose  phase  space  is 
considerably  reduced  compared  to  BaTiOs  and  KNb03.  Total  energy  calculations  of  the  AFD 
instability  showed  that  decreasing  the  volume  enhances  the  AFD  instability,  increasing  the 
corresponding  double-well  depth.  The  angle  characterizing  the  TiOe  octahedral  rotations  was 
found  to  be  about  a  factor  of  three  larger  than  observed,  however.  This  is  in  agreement  with  an 
earlier  LAPW  calculation,  and  after  extensive  tests,  we  ascribed  this  quantitative  discrepancy  to 
the  limitations  of  the  local  density  approximation.  [3] 

Even  though  no  further  phase  transition  is  observed  at  low  temperatures,  the  sharp  rise  in 
the  dielectric  constant  suggests  the  possibility  of  a  FE  instability  that  may  be  suppressed  by  the 
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zero-point  motion  of  the  atoms.  Such  an  instability  should  be  revealed  by  computing  the  full 
zone-center  dynamical  matrix  of  the  AFD  structure  (10  atoms  per  unit  cell).  We  focused  on 
investigating  the  influence  of  the  APD  transition  on  the  ferroelectric  instability  of  the  Fis  (with 
respect  to  the  cubic  structure)  ferroelectric  mode.  We  imposed  the  calculated  Fis  mode 
eigenvector  distortion  along  both  the  c-axis  and  a-axis  in  the  AFD  unit  cell  (the  TiOe  octahedral 
rotations  are  about  the  c-axis),  and  setting  the  twist  angle  of  the  TiOe  octahedra  at  the  1.4”  value 
measured  at  105  K.  We  found  that  the  AFD  structure  is  stable,  marginally,  against  the  imposed 
distortions.  [3]  Thus  large  ferroelectric-like  fluctuations  might  be  expected  in  the  AFD  structure, 
leading  to  the  increase  of  the  dielectric  function.  Since  the  AFD  structure  is  stable,  however,  with 
respect  to  the  ferroelectric  Fis  distortions,  the  crystal  never  actually  transforms  to  a  ferroelectric 
phase. 

C.  Impact  of  this  Work  on  Future  Research 

We  have  shown  that  while  the  AFD  structure  is  stable  against  subsequent  ferroelectric 
distortions,  it  is  only  marginally  so.  Thus  ferroelectric-like  fluctuations  might  be  expected  even  at 
low  temperatures,  resulting  in  the  observed  large  dielectric  function.  It  is  likely  that  future  work 
would  therefore  focus  on  these  fluctuations. 
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